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ABSTRACT 

We present new brightness and magnetic images of the weak-line T Tauri star 
V410 Tau, made using data from the NARVAL spectropolarimeter at Telescope 
Bernard Lyot (TBL). The brightness image shows a large polar spot and significant 
spot coverage at lower latitudes. The magnetic maps show a field that is predomi- 
nantly dipolar and non-axisymmetric with a strong azimuthal component. The field is 
50% poloidal and 50% toroidal, and there is very little differential rotation apparent 
from the magnetic images. 

A photometric monitoring campaign on this star has previously revealed V-band 
variability of up to 0.6 magnitudes but in 2009 the lightcurve is much flatter. The 
Doppler image presented here is consistent with this low variability. Calculating the 
flux predicted by the mapped spot distribution gives an peak-to-peak variability of 
0.04 magnitudes. The reduction in the amplitude of the lightcurve, compared with 
previous observations, appears to be related to a change in the distribution of the 
spots, rather than the number or area. 

This paper is the first from a Zeeman-Doppler imaging campaign being carried 
out on V410 Tau between 2009-2012 at TBL. During this time it is expected that the 
lightcurve will return to a high amplitude state, allowing us to ascertain whether the 
photometric changes are accompanied by a change in the magnetic field topology. 

Key words: stars: formation - stars: pre-main sequence - stars: rotation - stars: 
magnetic fields - stars: individual: V410 Tau 



1 INTRODUCTION 

T Tauri stars (TTS) are pre-main sequence late-type stars. 
As they contract towards the main sequence the increasing 
density in their interiors will lead to the development of a 
radiative core, if the star is sufficiently massive. A point of 
uncertainty concerns the evolution of their dynamos. In fully 
convective stars, or stars with small radiative cores, we ex- 
pect the dynamo to look fundamentally different to that of 
the Sun. We presume that the solar dynamo is concentrated 
in the boundary between the convective and radiative lay- 
ers, known as the tachocline. As the radiative core in the 
young star grows the dynamo should change from being dis- 
tributed throughout the convective zone to being confined 
to the tachocline. Observations of young stars with a range 
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of ages, masses and rotation rates will help to illuminate this 
process. 

V410 Tau is one of the most well-observed TTS. It is 
a weak-line T Tauri star (wTTS) which has already dissi- 
pated its disk, hence it does not display signatures of accre- 
tion characteristic of classical TTS (cTTS). I t has a vsini 
of 74 ± 3 km s _1 and a period of 1.872 days l|Stelzer et al.l 
2003), and is therefore an example of the young, rapidly ro- 
tating stars which are characterised by strong, large scale 
magnetic fields. Other , more evolved, ex amples include the 
10 Myr old TWA 6 liskellv et all 120081) 
Myr, iBarnes et al.ll200ll ' 
and AB Dor (50 Myr 



_ „, Speedy Mic (20 
lBarnesll2005l;lDunstone et alj|200( ' 
ICollier Cameron fc Donatil [200 



iDonati et ai1l2003al ; |jeffers et al.ll2007r i 

The strong fields in these stars have either been directly 
measured using Ze e man- Doppler imaging (ZDI. [s"eme]||l989l ; 
IDonati et al. 1 1 19971 . 1 2006), or their presence indicated indi- 
rectly e.g. by the large spot coverage on the surface. In the 
case of V410 Tau this field has revealed itself in the presence 
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of cool spots covering a large fraction of the photosphere. In 
fact, V410 Tau has exhibited one of the highest optical vari- 
abilities obs erved on a young s tar, with AV reach ing 0.6 
magnitudes i|Sterzer et alj|2003l ; IGrankin et al.ll2008l ). 

V410 Tau is the target of an ongoing photometric mon- 
itoring campaign which began i n 1981 llVrba et al. 1 19881 ; 
iHerbstl 1 19891 : IPetrov et all 1 19941 ; IGrankin et al.ll2008l ). For 
most of this time the lightcurve has been smooth, repeat- 
ing, with a clearly defined period and an amplitude between 
0.2 and 0.6 magnitudes. 

Two significant decreases in the amplitude of the 
lightcurve, where AV fell below 0.2, have been observed dur- 
ing this campaign. The first of these was between 1981-1983 
and the second began in 2005, reaching a minimum in 2007. 
Additional observations of V410 Tau dating from 1905-1987 
are available on photographic p lates in the ar c hives of the 
Sternberg Astronomy Institute. ISokoloff et al.l ([2008) have 
digitised this data and used it to investigate the variability 
over a longer period of time. This revealed a third decrease 
in lightcurve amplitude between 1963-1970, this time in the 
B-band. 

These changes in the lightcurves suggest a rearrange- 
ment of the cool spots. Given that starspots are a mag- 
netic phenomenon such a change is likely to be caused by a 
change in the magnetic topology. The fact that this change 
has been observed three times with intervals of approxi- 
mately 20 years may be evidence of a stellar cycle, similar 
to the 22-year solar cycle. The existence of magnetic cycles 
in active stars has been extensively investiga ted through 
chan g es in their starspo t coverage (see, e.g., iBerdvuginal 
120051 ; IStrassmeierl 120091 . and references therein). Possible 
short-period cy cles (of between 5-7 year s) have been foun d 
on V410 Tau bv lStelzer et al.l (120031 ) and lOlah et~aH l|2009h . 
using photometric dat a spanning severa l decades. 

Doppler imaging (|Vogt et al.l 1 19871 ; IStrassmeierl 120021 ) 
can provide more detail on the changes in spot distribu- 
tion. Maps of the spot dist ri butio n on V4 1 Tau have 
been made by | Joncour et alj (1 19941 ); iHatzesl (|l995al ) and 
iRice fc Strassmeiei ( 19961 ). The images have shown decen- 
tred polar spots and low latitude spots. Often the lower 
latitudes features are found to be confined with in particular 
longitude ranges (e.g. IRice fc Strassmeier|[l996l ). 

The flattening of the lightcurve suggests that a change 
in the distribution or number of photospheric spots has oc- 
curred, i.e., either there are fewer spots than previously 
or there is a more e ven longitudinal distribution of spots. 
IGrankin et alj (|2008l ) found a stronger correlation between 
the non-uniformity of the spot distribution and the ampli- 
tude of the lightcurve than the total spot area. This suggests 
that the reduction in lightcurve amplitude is related to a re- 
arrangement in the spot distribution, rather than a change 
in the number or area of the spots. 

To date, no magnetic images of V410 Tau have been 
published. This has motivated a long-term campaign, being 
carried out at Telescope Bernard Lyot (TBL), which aims 
to produce several spot and magnetic maps of V410 Tau 
between 2009-2012. We expect that during this time the 
lightcurve will return to the state observed between 1988- 
2004. In this paper we present the first set of images. Sec. 2 
discusses the evolutionary status of V410 Tau. Sec. 3 de- 
scribes the observations and data reduction; in Sec. 4 we 
describe the image production and present the brightness 



and magnetic images, and a new measurement of the dif- 
ferential rotation. Sec. 5 is a discussion of the implication 
of the results for the dynamo theory in pre-main sequence 
stars. 



2 PHYSICAL PARAMETERS OF V410 TAU 

The evolutionary state of V410 Tau is somewhat uncertain. 
The v sin i and period quoted above indicate a minimum 
stellar radius of 2.8 R©. Widely varying effective tempera- 
tures have been reported in the literature, incl uding 4700 K 
(e.g.. iBertout et al.l [2007h and 4060 K (e.g., ISestito et ail 
120081 ). In this work we adopt a temperature of 4500 ± 100 K, 
in better agre e ment with the B - V value of 1.16 from 



Grankin et al 



Bessell et alj 



(2008) and the temperature calibrations of 



1998). 



With this temperature and minimum radius we 
can estimate th e min imum luminosity t o be 3.4 Lq . 
IWichmann et al.l (|l998l ) and IBertout et al] (|l999l ) g ive a 
range of distances of 116-157 pc, based on Hipparcos l|ESAl 
Il997h measurem ents. Assuming the lowest observed V mag- 
nitude (10.52) in lGrankinl(|l999l ) is the unspotted V magni- 
tude, and a distance of 157 pc gives a luminosity of 2.4 Lq 
(if A v = 0) , well below the minimum luminosity calculated 
above. 

The discrepancy suggests that even when we observe 
the star at its brightest, it still has significant spot coverage. 
In other words, the true "unspotted" mag nitude is brighte r 
than 10.52. This possibility was explored in lGrankm] (l999), 
who calculated an absolute unspotted V mag of 10.236, 0.28 
mag brighter than their minimum observed magnitude. This 
would mean that the star has an even distribution of small 
spots across its surface at all times that are not detected 
by Doppler imaging, as they are smaller than the resolution 
limit. 

Using an inclination of 70 ± 10° (obtained in the imag- 
ing process, see Sec. I4.2J1 gives a radius of 3.0 R0, and with 
a temperature o f 4500 K indicates (using the models of 
ISiess et al.lboOOh that V410 Tau has a mass of 1.4 ± 0.2 M e 
and an age of 1.2 ±0.3 Myr. This is much you nger than pre- 
viously indicated, e.g. in lBertout et al.l (|2007l ). where an age 
of 6 Myr was given for V410 Tau. The parameters given in 
IBertout et al.l (|2007-) are consistent with a radius of 1.8 R0, 
which is incompatible with the vsini and period. 

These param eters indicate (again using the models of 
ISiess et al.l l2000l ) that the radiative core has a radius of 
0.1R*. The uncertainties in the temperature and inclina- 
tion allow the radiative core to have a radius between 
0.0— 0.28R* , meaning that the star could be fully convective. 

Some interesting questions arise when we consider the 
possibility of a stellar cycle, particularly as we expect that 
a star which is fully convective or has a small radiativ e core 
has close to 'solid body' rotation (jBarnes et al.| [2005). 

Differential rotation is crucial to the solar magnetic 
cycle, as it allows the development of toroidal field from 
poloidal field by shearing the field lines (Parker 1955). This 
may not rule out the possibility of cyclical dynamos exist- 
ing wi thout differential rotation. For example. iRiidiger et al.l 
(2003) modelled a 2 dynamos and found that oscillatory so- 
lutions appear in turbulent spherical shells in exceptional 
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cases. However, they conclude that cyclical activity is al- 
most always indicative of differential rotation. 

Previous attempts to detect diffe rential rotation on 
V410 Tau have found it to be very small. iRice fc Strassmeier] 
obtained a value for the differential rotation param- 
eter An of 0.003, around 1/20 of the solar value. If the star 
does indeed have a cycle it appears this must arise with very 
little differential rotation. 



3 OBSERVATIONS AND DATA ANALYSIS 

Spectroscopic observations of V410 Tau were made at the 
2-m Telescope Bernard Lyot (TBL) on Pic du Midi in the 
French Pyrenees between 02-17 January 2009. During this 
time, 40 spectra were taken using the spectropolarimeter 
NARVAL (the twin of ESPADONS on the Canada France 
Hawaii Telescope), with four sub-exposures taken in dif- 
ferent polarimeter configurations each time. This allows 
spurious polarisatio n signals to be removed to first order 
l|Donati et al.|[l997l ). The integration times (for the individ- 
ual sub-exposures) were 600s, except for the observations 
taken on 10 January 2009 which had integration times of 
800s. 

The full wavelength range of the spectra was 370 — 
1000 nm, taken with a spectral resolution of 60000 and a 
velocity pixel step of 2.6 km s _ , reduced to approximately 
1.8 km s _1 as the spectrograph slit was tilted with respect 
to the CCD pixels. The full observing log is given in Tab. [1] 
In the table the rotational phase of V410 Tau at each obser- 
vation is given with respect to the ephemeris, 



JD(Hel.) = 2454832.58033 + 1.87197R, 



(1) 



where R is an integer number o f rotations. This us es the 
same period and zero phase as IStelzer et afl (|2003l) . The 
resulting phase coverage is good, covering around 75% of 
the surface and there is close phase correspondence between 
several of the observations providing a good opportunity to 
measure differential rotation (see Sec. I4.4|l . 

The data was reduced using Libre ESpRIT, a dedi- 
cated pipeline for the reduction of spectropolarimetric data 
taken at TBL and CFHT. This package automatically ex- 
tr acts the unpo l arised and polarised spectra, as described 
m iDonati et all (|l997h . To improve the effective signal-to- 
noise ratio of the data least-squares deconvolution (LSD, 
IDonati et ai1ll997T l was applied to each spectrum, i.e., a line 
list was deconvolved from each spectrum, revealing in each 
case the underlying rotationally broadened profile. Mean in- 
tensity (Stokes I) and circularly polarised (Stokes V) profiles 
were produced in this way. The line list wa s created using an 
ATLAS9 model atmosphere l|Kuruczlll99l L using a spectral 
type of K5, and contained just under 9000 absorption lines 
in total. The Stokes I profiles show distortions characteristic 
of photospheric spots and Zeeman signatures were detected 
in most spectra, with a typical peak-to-peak amplitude of 
0.2%, i.e. 10 times the noise level. Example Stokes I and V 
profiles are shown in Fig. Q] 

We used the I and V profiles to directly calculate the 
magnetic field along the line-of-sight. This field is calculated 



using 



B«(G) 



Table 1. Spectra taken at NARVAL/TBL. The dates are all in 
January 2009. The observing times shown indicate the beginning 
of the integration. Column 4 shows the maximum signal-to-noise 
ratio per CCD pixel. Column 5 is the rms noise level of the po- 
larised LSD profiles. 
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2.12 
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2.10 
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02 


34.38435 


21:08:39 
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2.08 


0.9637 


02 


34.41572 


21:53:50 
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2.07 
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02 
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158 
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03 
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19:56:30 
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2.31 
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04 
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2.17 
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21:08:28 
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10 
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5.2286 
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10 
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2.31 
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11 


43.29158 
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135 
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11 
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141 
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5.7559 


11 


43.38575 
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46.36540 
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145 


2.48 
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14 
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2.49 
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3.37 
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49.26621 


18:20:13 
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3.01 
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17 
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19:05:25 


102 


4.03 
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17 


49.32898 


19:50:37 


120 


3.16 
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17 


49.36037 


20:35:49 


139 
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17 


49.39175 


21:21:01 


136 


2.71 


8.9813 



-2.14 x 10 1 



/ vV(v)dv 



AogC / [Ic - I(v)]dv 



(2) 



from IDonati et all l|l997t) where v is the velocity in the star's 
rest frame, Ao the central wavelength of the LSD profile in 
nanometres, c is the velocity of light in vacuum, g e g is the 
value of the mean Lande factor of the LSD line and I c is the 
continuum level. This is plotted in Fig. [2] Near phase 0.0 
the field peaks at a value of around 200G. The value falls to 
a plateau between phases 0.25 — 0.5 where the field is about 
150 G in the other direction. Although the individual error 
bars are large the plot reveals the rotational modulation 
of the longitudinal field, as the pattern of variability just 
described is repeated throughout the eight rotations that 
were observed. 

Photometric observations of V410 Tau were obtained 
with two different photoelectric telescopes at different lo- 
cations in Crimea, Ukraine: the 0.6-m T-60 on Mt. Koshka 
equipped with UBV(R)c filters for the Johnson-Cousins sys- 
tem, and the 1.25-m AZT-11 in Nauchny equipped with the 
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-100 100 

Velocity (km/s) 

Figure 1. Example Stokes I (bottom) and V (top) profiles at 
phase 7.331. The V profile has been multiplied by fifteen and 
shifted upward by 1.03. 
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Figure 2. Line-of-sight magnetic field calculated using Eq.[2] The 
black triangles represent rotations — 5 (as listed Tab. [TJ while 
the red crosses represent rotations 6 — 8. 

similar filters. Altogether, 41 UBV(R)c data points were col- 
lected for this star between 03 August 2008 and 26 January 
2009. All measurements were made differentially with re- 
spect to a nearby comparison star and a check star in the 
the following sequence: CK-S-C-V-C-V-C-V-C-S-CK, where 
CK is the check star, C the comparison star, S the sky back- 



ground, usually between the comparison and the variable 
star, and and V is the variable itself. Integration time was 
usually set to 10s for all filters. The data reduction was 
based on nightly extinction coefficients obtained from a set 
of standard stars. The rms error of a single measurement 
in the instrumental system is about 0.005 mag in BVR and 
0.01 mag in U. Only the V-band lightcurve is used in this 
work. 



4 DOPPLER AND MAGNETIC IMAGES 

4.1 Zeeman-Doppler imaging process 

Doppler imaging produces maps of cool photospheric spots 
by tracking the distortions that the spots cause to rotation- 
ally broadened spectral lines as the star rotates. Zeeman 
Doppler imaging extends this technique to mapping mag- 
netic fields by applying a similar process to the circular po- 
larisation profiles. 

Details of the Zeeman Doppler i maging process have 
bee n given in several p apers including iDonati et all (|l997t ) 
and lDonati et all (|2006T ) . The code finds the best-fitting im- 
ages using x 2 minimisation using the I and V profiles as 
constraints. This cannot provide a unique solution, as the 
spatial resolution of the image is always limited to some ex- 
tent. For this reason, entropy maximisation is used as an 
additional constraint, i.e. the information content of the im- 
age is minimised. 

The ZDI process usually begins with a spotless star and 
a null magnetic field. At each step in the iteration the I 
and V profiles corresponding to the current image are cal- 
culated. The local line profiles are described using Gaussian 
profiles, which has been found to be adequate in the case of 
rapidly rotating stars where the shape of the local line pro- 
file is less important. The central wavelength of the mean 
line is 620 nm, the Lande g factor of 1.2 and a FWHM of 
9.7 km s~ . A linear limb-darkening law with a coefficient 
of 0.75 was used. The code explores the parameter space, 
finding stellar configuration which gives Xr ~ 1-0 while si- 
multaneously maximising the image entropy. 

4.2 Brightness Image 

The free parameters in the production of the brightness im- 
age are the spot-filling factors of each pixel of the image 
(Colli er" Cameronll 19921 ). i.e., the fraction of each pixel which 
contains spot. 

V410 Tau is a well-studied star, therefore the stellar 
parameters, including the usini, period, radial velocity, in- 
clination and differ ential rotation, have been measured in 
previous work (e.g . jHatz es 1995b; Rice fc Strassmeierlll996l ; 
IStelzer et all 120031 ) . We can confirm these values using v 
minimisation before producing the final images. To fit the 
stellar parameters 1000 image pixels were used, while the 
final images were produced using 10000 pixels. The final 
reconstruction was carried out using a radial velocity of 
16.8 km s _1 , projected rotational velocity v sin i — 75 km s" 1 
and inclination i = 70° (consistent with the values in the in- 
troduction). The fitting of the differential rotation and pe- 
riod is described in Sec. 14.41 

Fig. [3] shows the measured intensity profiles and their 
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Figure 4. Doppler image of V410 Tau shown as a flattened polar 
projection. The colour scale represents the spot-filling factor. The 
circles represent lines of latitude separated by 30°, the largest 
circle is at latitude —60° and the bold line is the equator. The 
numbers around the outside are rotational phases and the outer 
tickmarks show phases at which observations were made. 

fits, while the Doppler image is shown in Fig. [4] It reveals a 
polar spot and groups of spots at lower latitudes. The total 
spot filling factor on the visible regions of the star is 8%. 
The lower latitude spots are concentrated in the latitude 
band between — 30° and are fairly evenly spread longitu- 
dinally. The relationship between the spot distribution and 
the optical variability will be discussed in Sec. 15.11 

4.3 Magnetic images 

The poloidal and toroidal components of the magnetic field 
are described using a spherical harmonic expansion with 
each component of the field (the radial, non-radial poloidal 
and toroidal field components respectively) expressed in 
spherical c oordinates. The equa tions describing the field are 
set out in iDonati et al.1 l|2006ft . The advantage of using a 
spherical harmonic decomposition, rather than allowing each 
pixel to vary freely, is that the spherical harmonic descrip- 
tion gives rise to more physically realistic solutions. For ex- 
ample, the code will straightforwardly return a dipole, which 
is unlikely when the field values in individual pixels are left 
as independent variables. 

As the spatial resolution depends on the ratio of the 
width of the absorption lines to the velocity resolution in 
the spectra, the high v sin i of V410 Tau translates to a high 
resolution on the stellar surface. This allows the field to be 
fitted up to I = 15. Fig.[5]shows the resultant V profiles and 
the fits. 

The radial, azimuthal and meridional components of the 
field are shown in Fig. [6] These images show a strong, ap- 
parently complicated magnetic field, with individual regions 
having strengths of up to 1 kG. The average field integrated 
over the surface of the star is 490 G. The radial component is 
complex with field covering a large proportion of the surface, 



consisting of many small regions of magnetic field and sev- 
eral switches in polarity. The azimuthal field covers almost 
the entire surface of the star, with field strengths between 
300 — 1000 G covering more than half the visible surface. 
There is a monolithic region of positive polarity completely 
covering one hemisphere, with two smaller regions of neg- 
ative polarity in the opposite hemisphere. The meridional 
field is weaker than the radial and azimuthal components 
and shows positive and negative regions near the pole. 

We can quantify the complexity of the field by calcu- 
lating the fractional strengths of different components. The 
field is evenly split between the poloidal (49.2%) and toroidal 
(50.8%) components. The poloidal is dominated by the dipo- 
lar component, accounting for 60% of the total field. The 
poloidal field is nearly completely non-axisymmetric, with 
5.3% of the poloidal field being axisymmetric. The implica- 
tions of the magnetic topology in relation to the fields in 
young stars in general will be discussed in Sec. 15.21 

4.4 Differential Rotation 

We carried out a measurement of differential rotation on 
V410 Tau by repeating the image production process de- 
scribed above, while incorporating a solar differential rota- 
tion law of the form 

n(0) = Q cq - Afi sin 2 e, (3) 

where £1(8) is the angular velocity as a function of latitude 
(8), Sl cq is the angular velocity at the equator and AQ is the 
difference between the angular velocities at the equator and 
the poles. 

To constrain the values of Q eq and AQ, we varied their 
values while holding the information content of the images 
constant, and produced a 2D plot of Xr- F° r the Stokes I 
profiles the percentage spot coverage was fixed and for the 
Stokes V profiles the the magnetic field strength (B a i m ) was 
held constant. When the I profiles were used we were unable 
recover the differential rotation parameters. The difficulty 
here arises because the spot distribution is essentially a polar 
spot and a band of lower latitude spots between — 30° , min- 
imising the differences between the resultant images when 
differential rotation parameters are introduced. 

However when the V profiles alone were used in the 
minimisation it was possible to fit a paraboloid to the \ 2 
plot. The minimum of this paraboloid identifies the Q eq and 
AQ favoured by the data. When the code was allowed to 
converge with a fixed field of 400 G we obtained the \ 2 
surface shown in Fig. [7] The minimum of the paraboloid 
fitted to this surface lies at Q cq — 3.357 ± 0.001 rad day -1 
and AQ = -0.007T0.006 rad day" 1 . This is compatible with 
solid body rotation as photospheric shear is less than 30% 
the solar value at the 3-sigma level. When different values 
of B a i m are used we recover very similar values of Q cq , AQ, 
and their uncertainties. 

4.5 Balmer Line Analysis 

The Balmer line emission provides another probe of the 
surface activity of TTS. The Ha profile in wTTS is con- 
sidered to be a superposition of a photospheric absorption 
component and a emission component originating in the 
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Figure 3. The I profiles produced using LSD are plotted in black and the fits are plotted with a thin red line. The number of rotational 
profiles completed since the ephemeris is given on the right. 
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Figure 6. The magnetic images of V410 Tau, displayed using 
the same projection as Fig. [4] The colour scale represents the 
magnetic field in units of Gauss. 



chromosphere and corona (cTTS h ave additional contribu - 
tions from the disc and jet), see e.g Kurosawa et al.l (2006); 
iBarrado v Navascues fe Martini (|2003l ). The Ha profiles in 
V410 Tau, taken from the spectra listed in Tab.[T] show con- 
siderable variability, with the line appearing in both absorp- 
tion and emission. Example profiles showing this variation 
are given in Fig. [8] one showing the line appearing in emis- 
sion, the other when it is in absorption. In order to highlight 
the phase- dependence of the variability Fig.|5]shows the dy- 
namic spectrum. This is created by stacking the profiles in 
order of phase and displaying as a colour plot. 

The dynamic spectrum shows emission and absorption 
features. Ha is most strongly in emission in the phase range 
0.1-0. 5, suggesting the prese nce of chromospheric active re- 
gions (jMekkaden et alj|2005l ) on that hemisphere. These ac- 
tive regions may be connected to the large region of az- 
imuthal field on that hemisphere (see Fig. [6} . 

This emission feature is crossed by an absorption tran- 
sient that is centred around phase 0.4. There is a similar 
absorption transient between phases 0.9-1.05. When such 
absorption transients have been observed in other stars they 
have be en identified with prominences crossing the stel- 
lar di sc l|Collier Cameron fe Robins on 1989; Dunstone et al.l 
2006). The absorption transients appear to be passing 
through the spectral line more quickly than the spots (i.e. 
their slopes in Fig. [9] are low), suggesting that they are 
not ver y close to the surf ace. Calculating the variance 
profile ( Johns fe Basril Il995l ) reveals a symmetrical profile 
(not shown), with a full- width half- maximum (FWHM) of 
130 km s _1 , and shows variance out to ±100 km s _1 . 

The small amount of variability at velocities larger than 
v sin i appears to suggest that a large fraction of the vari- 
ability is due to chro mospheric activity. Whi le this has been 
observed previously (|Fernandez et all 120041 ') , we currently 
favour higher-lying features due to the flat slopes of the 
absorption transients in Fig. [9] The absorption feature be- 
tween phases 0.90 and 1.05 was observed throughout the 
observation run, i.e. for 8 stellar rotations, suggesting the 
presence of a stable feature in corotation with the star, 
such as a prominence. If this is the case we can estimate 
the height of the prominences from the slopes of the fea- 
tures in Fig. [9] The absorption transients cross the stel- 
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Figure 8. Example Ha profiles, showing the line in emission at 
phases 1.522 (top) and in absorption at phase 0.998 (bottom). 



lar disc in roughly 0.15 of a rotation. Assuming the promi- 
nences are in co-rotation with the star this corresponds to 
a distance of 2.5 R, cos 9 from the rotation axis, where 6 is 
the stellar latitude. The co-rotation radius is 2.4 R*, sug- 
gesting they may lie above co-ro tation. The equations in 
Ijardine fc van Ballegooiienl (|200 5) allow stable prominences 
up to a height of 7.4 R„. 



5 DISCUSSION 

5.1 Optical Variability of V410 Tau 

At the present time V410 Tau is displaying a low ampli 
tude optical lightcurve, which we can atte mpt to relate 



to the image in Fig. [4] Previous images (IJoncour et al 



19941 ; Hatzcs 1995a; Ricc fc Strassmeiedl 19961 ; ISchmidt et al 



20051 ') showed decentred polar spots and low latitude spots 
confined to particular longitude ranges. In Fig. [4] the centre 
of the high latitude spot is very close to the pole and the 
lower latitude spots are found at all longitudes, leading to 
less variability as the star rotates. 

This variability can be quantified more precisely by 
comparing with the photometric data and calculating the 
flux differences that result from the spot distribution in 
Fig. 2] The photometric measurements in the V-band made 
at the Crimean Observatory and described in Sec. [3] are 
shown in Fig. 1101 A synthetic lightcurve has been overplot- 
ted. This was produced by integrating the flux over the vis- 
ible surface of the star at each phase. The spot-filling fac- 
tors were converted to temp eratures, using Tphot — T S pot — 
1400 K (|Petrov et all Il994h . and then into fluxes, which 
are then averaged over the visible surface of the star at 
each phase (taking limb-darkening into account) and finally 
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Figure 9. Dynamic spectrum of the Ho profiles. Red represents 
emission and blue is absorption and the normalised intensity 
range is between 0.7 and 1.3. The vertical lines mark the po- 
sition of the i; sin i (75 km s — 1 ). The phase range on the vertical 
axis is between 0.1-1.1 to better display the continuation of the 
absorption feature around phase 1.0. 



converted to a AV. This gives a smoothly varying func- 
tion with a peak brightness at phase 0.85 and a minimum 
at phase 0.25. The peak-to-peak amplitude of this curve is 
0.038 which is roughly consistent with the photometric mea- 
surements. 

We believe that the outlying point at phase is asso- 
ciated with a flare. The contemporaneous U-B and B-V 
colours provide additional evidence for this event. Discount- 
ing this point, the dispersion in the measured data points 
at this time is about 0.02. The mean V magnitude mea- 
su red during this epoch (10.85) is similar to that measured 
in iGrankin et al.l ()2008h where < V> lies between 10.8 and 
10.95, further supporting the assertion that the change in the 
lightcurve is associated with a rearrangement in the spots. 



5.2 Magnetic Nature of V410 Tau 

Zeeman Dopplcr images have previously been created 
for four other TTS, all of which are cTTS. These are 
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Figure 10. V magnitude as a function of rotational phase. The 
circles show data taken at the Crimean Observatory. The aster- 
isks show the lightcurve calculated by integrating the flux at each 
phase from the image shown in Fig. [4] The latter has been shifted 
to have the same mean V magnitude as the photometric measure- 
ments. 



BP Tau (jDonati et alj l2008h. V2129 Oph (IDonati et all 
12007t h CR Cha and CV Cha l|Hussain et all2009h . There is a 
trend observed in M dwarfs, that fully convective stars have 
large proportions of their fields concentrated in the lowest I 
and m modes and stars with radiative cores hav e more power 
concentrated in the hig her I and m modes l|Morin et al.l 
l200Sl ; lDonati et alJliooa j. BP Tau is a fully convective cTTS 
with a field which is strongly axisymmetric and dipolar, sim- 
ilar to fully convective main sequence stars. V2129 Oph, 
which has a radiative core, is dominated by an octupolar 
component which is approximately three times larger than 
the dipole. CR Cha and CV Cha are also partially convec- 
tive and have significant power in higher I and m modes. 
The field of V410 Tau has a strong dipolar component, in 
common with fully convective stars, despite the fact that 
it appears to be partially convective (see Sec.[2|. The large 
fraction of the field confined into the dipolar modes is signif- 
icant, given that the high v sin i of V410 Tau means that we 
are sensitive to the higher orders in the spherical harmonic 
decomposition. 

Fig. [TTJ compares the magnetic topology of V410 Tau 
with BP Tau and V212 9 Oph, using the conventions in 
iDonati et ail ([2008) and iMorin et al.l (|2008l ). who investi- 
gated the magnetic topology of M dwarfs. While we cannot 
infer trends from this plot from only three data points it 
allows direct comparison between these TTS. V2129 Oph 
and V410 Tau have similar physical structures, with closely 
matched masses and similar radiative core sizes (~ 0.2 R„ 
for V2129 Oph). However, there are significant differences 
between the fields of the two stars. The field of V410 Tau 
is more strongly poloidal, and is much less axisymmetric. 
ft is rotating more rapidly than the cTTS BP Tau and 
V2129 Oph (due to the absence of disk-braking), but this 
does not lead to commensurately larger fields as it is the 
fully convective BP Tau which generates the largest mag- 
netic field, despite its lower rotation rate. It appears that 
the depth of the convection zone plays a role in determin- 
ing the field strength, although we do appear to observe the 
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Figure 11. Plot of the mass of TTS against rotational period. 
The size, colour and shape of symbols represent the field strength, 
configuration and axisymmetry of the fields respectively. Larger 
symbols indicate a larger field, red stands for 100% poloidal and 
blue stands for 100% poloidal. A sharp five-pointed star represents 
a completely non-axisymmetric poloidal field and a decagon rep- 
resents a completely axisymmetric poloidal field. The dotted line 
shows the approxima te position of the fully convective limit for a 
1 Myr star using the lSiess et alj l|200Ct ) models. 



sharp transition in field topology that is seen in M dwarfs. 
We anticipate adding more stars as the data on TTS become 
available. 

Qualitatively, the field of V410 Tau most resem- 
bles older, rapidly rotating stars such as AB Dor 
dDonati fe Collier Cameronll 19971 ; IDonati et alll 19991 . l2003bl ; 
iHussain et al.ll2007l ) , rather than the classical TTS discussed 
above. AB Dor is a 50 Myr old (i.e. post-T Tauri) K0 star 
with a usini of 90 km s . The similarities between the 
magnetic maps of the two stars are monolithic regions of 
azimuthal field, almost forming rings around the pole and 
complex radial fields with many switches in polarity. 

AB Dor and V410 Tau are broadly similar in the sense 
that they are young and rapidly rotating. However V410 Tau 
has almost solid body rotation while AB Dor has significant 
differential rotation, with AQ (s ee Eq. having values be- 
tween 0.043 - 0.074 rad d ay" 1 (ICollier Cameron fc Donatil 
120021 ; IDonati et al1l2003l lal; Ijeffers et al.ll2007l ). 

Usually, strong azimuthal field is associated with dif- 
ferential rotation. For instance an aQ dynamo gives rise 
to toroidal field through differential rotation shearing the 
poloidal field. The field topology of V410 Tau is therefore 
somewhat difficult to understand. 

While the qualitative similarity of the field with those 
of partially-convective stars such as AB Dor suggests that 
V410 Tau has a radiative core, we cannot rule out the possi- 
bility that V410 Tau is fully convective, as Sec. [2]suggested. 
In this case the dipolar component of the field follows the 
trend observed in other TTS and M dwarfs. However, the 
question of the origin of the azimuthal component remains, 
due to the lack of observed differential rotation. 

Further magnetic and brightness maps will follow from 
this campaign on V410 Tau. We also have spectra of several 
other TTS which will be used to produce magnetic images. 
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Using these data, we expect to gain a better insight into the 
form and evolution of dynamos in pre-main sequence stars. 



6 CONCLUSIONS 

The Doppler image of the weak-line T Tauri star V410 Tau 
in January 2009 shows a large polar spot and a several 
lower latitude spots spread around the latitude band be- 
tween — 30°. This spot distribution can explain the rel- 
atively flat V-band lightcurve currently being observed on 
the star as, unlike during previous Doppler imaging cam- 
paigns, the spots are not confined to a particular longitude 
range. The polar spot on the image is centred on the pole, 
unlike on previous images where it was shifted from the pole. 

The Ha profile is complex and highly variable. Fast- 
moving absorption transients in the line may indicate the 
presence of high-lying circumstellar material, such as sling- 
shot prominences. 

The magnetic field is predominately dipolar and non- 
axisymmetric. It has a complex radial field as well as a strong 
azimuthal field with large monolithic regions. The field re- 
sembles zero-age main sequence stars such as AB Dor, but 
unlike that star it has very weak differential rotation. The 
rotation rate and depth of convection zone appear to be the 
most important parameters in determining the strength and 
configuration of the magnetic field 

Future Zeeman Doppler images planned as part of a 
campaign on Telescope Bernard Lyot (TBL) will reveal 
whether changes in the amplitude of the optical lightcurve 
are accompanied by significant changes in the spot distribu- 
tion or magnetic field topology. 
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